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The kinematic regime of the magnetic dynamo neglects the backreaction of the 
magnetic field on the flow. For small magnetic diffusivity, in the early stage of 
evolution, there is an ideal phase where dissipative effects can also be neglected. We 
estimate the magnitude of the energy dissipation term (Ohmic heating), taking into 
account differential constraints on chaotic flows. We find that the period of ideal 
evolution is roughly doubled over an estimate without constraints. The helicity 
generation terms are exponentially smaller than the energy dissipation, so that large 
quantities of energy are dissipated before any helicity can be created. Helicity flow 
is exponentially larger than net helicity generation. The constraints also lead to 
the existence of a singular initial condition for the magnetic field for which sizable 
amounts of helicity can potentially be created. 
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I. INTRODUCTION 

For many three-dimensional flows an embedded magnetic fleld will grow exponentially; 
this phenomenon is referred to as a dynamo [1]. The kinematic dynamo problem consists 
of studying the induction equation for the fleld on the assumption that it does not react 
back on the flow (the Lorentz force is neglected). This is justifled as long as the fleld is 
small enough compared to inertial forces. A good estimate of the domain of validity of the 
kinematic assumption is thus of paramount importance. 

It was recently pointed out by Schekochihin et al. [2] that the details of how the nonlinear 
Lorentz force term appears in the equations of motion are important. The Lorentz tension 
force only involves the gradient of the magnetic fleld along itself (the parallel gradient), and 
in smooth chaotic flows or turbulent flows at large Prandtl number that gradient is much 
smaller than its perpendicular counterpart. Thus, the onset of the reaction force is delayed 
as compared with an estimate based on the scale and magnitude of the magnetic fleld. 

In this paper we investigate the domain of validity of the ideal evolution assumption, 
that the fleld is passively advected by the plasma with resistive effects being negligible. 
Because the resistivity (magnetic diffusivity) is assumed very small in dimensionless terms 
{i.e., the magnetic Reynolds number is very large, typically 10^-10^^; for example, it is of 
order 10^^ in the solar corona), there is a sizable period of evolution where the scale of 
the magnetic fleld is such that Ohmic dissipation is unimportant. In a flow with chaotic 
Lagrangian trajectories the exponential stretching of fluid elements leads to growth of the 
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amplitude of the magnetic field, so that the magnetic energy grows very rapidly during the 
period of ideal evolution. But in incompressible flows there is a corresponding exponential 
decrease via folding of the scale of magnetic field variations. These small scale variations in 
the field cause the Ohmic dissipation term to become large, and the evolution is no longer 
ideal. A convenient method of estimating the importance of the nonideal (resistive) term is 
to compare the rate of growth of magnetic energy to that of the power dissipated by Ohmic 
heating. This was done in Refs. [3, 4] where it was found that, because the power dissipation 
involves extra gradients of the magnetic field, the ratio of energy over power (the dissipation 
time scale) decreases exponentially in time at a rate — 4A, where A is a typical value of the 
Lyapunov exponent of the flow — the mean rate of exponential stretching of fluid elements. 

In the present paper we will show that these estimates are overly pessimistic. Because 
of the existence of differential constraints on the characteristic directions and amplitudes of 
stretching, the dissipation time scale decreases exponentially at a typical rate —2 A, implying 
that the period of ideal evolution is twice as long as previously thought. These differential 
constraints originated in the work of Tang and Boozer[5] on the advcction-diffusion equation , 
and Thiffeault and Boozer [6] found a new constraint valid in three dimensions that we 
apply here to the kinematic dynamo. Methods from differential geometry were used to show 
that the derivatives along the characteristic directions of expansion of the stretching rates 
cannot be completely independent and have to satisfy constraints. The theory of differential 
constraints was recently expanded and put on a more rigorous footing by Thiffeault [7]. The 
constraints have been tested numerically to a high degree of precision for various maps and 
flows [5-7]. 

After reviewing aspects of stretching in chaotic flows in Section II, we calculate the 
energy dissipation time scale in the manner of Boozer [3] in Section III and introduce the 
necessary notation along the way. The effect of differential constraints on the dissipation time 
scale is derived in Section IV. In Section V we investigate an intriguing consequence of the 
differential constraint, namely that the exists a singular initial condition for which the energy 
dissipation time scale does not decrease exponentially; we speculate on the physical meaning 
of such an initial condition. Section VI is devoted to deriving estimates of the amount of 
helicity created during the period of ideal growth, both with and without constraints and 
for the singular initial condition. Finally, in Section VII we offer some general comments on 
the results of the paper and discuss future work. 



In this section we review the properties of stretching in chaotic ffows, as embodied by a 
metric tensor in Lagrangian coordinates that measures the total deformation experienced 
by a fluid element as it is advected and stretched by the flow. 

For a given smooth velocity held v{x, t), the trajectory of a fluid element is the solution 

to 



where x{a, t) is the position at time t of a fluid element that started at a at t = 0. The initial 
condition labels a are called Lagrangian coordinates, and the x are the Eulerian coordinates. 
(The Eulerian coordinates are "flxed" in space and can be regarded as ordinary Cartesian 
coordinates.) The function x = x{a,t) is thus the transformation from Lagrangian (a) to 
Eulerian {x) coordinates. For a chaotic flow, this transformation gets extremely complicated 



II. STRETCHING IN CHAOTIC FLOWS 
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as time evolves. 

The Jacobian matrix of the transformation x{a,t) is M'g := dx^/da'^. The Jacobian ma- 
trix is a precise record of how an initially spherical fluid element is rotated and stretched 
into an ellipsoid by the flow. We are interested in describing stretching and its direc- 
tion in the reference frame of a fluid element (Lagrangian frame), not the absolute ro- 
tation in Eulcrian space, so we construct the metric tensor or Cauchy-Green strain ten- 
sor [8], Qpq := Ylt=i M^pM^qj which contains only the information on the stretching of fluid 
elements (this can be seen, for example, by performing a singular value decomposition 
of M [7, 9, 10]). The metric tensor is a symmetric, positive-definite matrix, so it can be 
diagonalized with orthogonal eigenvectors {u, m, s}, and corresponding real, positive eigen- 
values {A^, A^, Ag}; we can thus write 

The A's and {u, m, s} are functions of a and t (but see below). The A's are called coefficients 
of expansion. Without loss of generality, we order them such that A^ > Am > Ag (we assume 
nondegeneracy of the eigenvalues) . The label "u" indicates an unstable direction: after some 
time, we have Au ^ 1, growing exponentially for large time. Such a direction exists if the 
flow has chaotic trajectories, which we assume to be the case. The label "s" indicates a 
stable direction: after some time, we have Ag <^ 1, decreasing exponentially for long times. 
The intermediate direction, denoted by "m" , can grow or decrease exponentially at a rate 
somewhere between that of Au and Ag; for steady flows {i.e., when v is independent of t), Am 
has only algebraic (nonexponential) behavior [11]. The determinant of the metric is denoted 
by g; if v is incompressible (V • v = 0), we have ^/g = Au AmAg ~ 1. 

The eigenvectors {u, m, s} are called characteristic directions of stretching. It is well- 
known that these characteristic directions converge exponentially in to time-asymptotic 
vahics that depend only on the initial condition a (and the initial time for unsteady flows) [7, 
12]. We assume that these directions (on the trajectory of interest) have converged. It is 
also well-known that the Lyapunov exponents converge far slower than the characteristic 
direction. We do not assume that the finite-time Lyapunov exponents Xfj{a,t), defined by 

A^:=^logA^, (3) 

have converged or even that their limit as f — > oo exists; we merely suppose that t is 
sufficiently large to have Au ^ 1 and Ag <^ 1. This difference in magnitude between the 
unstable and stable directions is the basis for the approximations used in this paper. 

III. ENERGY AND CURRENT ALONG TRAJECTORIES 

Equipped with the metric tensor in diagonal form (2), we can now compute the evolution 
of physical quantities, such as the energy and current, along fluid trajectories. This was 

investigated in Rcfs. [3, 4], where the relative growth rates of energy, parallel current, etc. 
were derived. In the present section we revisit these results; in Section IV we will show that 
they must be modified to account for differential constraints. 

The evolution of a magnetic field in resistive MHD is governed by the induction equation. 



= V X (v X B) + —V^B 



(4) 
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where B is the magnetic field, rj is the resistivity, and //q is the permeabihty of free space. 
With the help of the chain rule and the metric tensor, we can transform the magnetic 
induction equation (4) from Eulerian coordinates, x, to Lagrangian coordinates, a. 



d_ 
di 



na,t) = Y: ' ^ 



p,<i- 



Ho daP 



(5) 



where V" / ^ '■= Ylii i^'^Yi^^ the magnetic field in the Lagrangian frame, and gi^^ := 
[g-'^)Pi. The la subscript on d/dt is a reminder that the time derivative is taken holding a 
constant, as opposed to the derivative in (4) which has x fixed. 

Equation (5) is simply a difi^usion equation with anisotropic difi^usivity tensor {r] / fio) g^'^ . 
By construction, the velocity v has dropped out of the equation entirely. When rj = 0, 
we have the well-known Cauchy solution b = 6(a), independent of time, so that b is the 
magnetic field &t t — 0. This is the classic result that in ideal MHD the magnetic field is 
frozen into the plasma. 

In the absence of diffusivity, the magnetic field is amplified in a chaotic fiow, as can be 
seen by considering the magnetic energy Eb '■= (l/2/io)-B^ = g^qlP If /2iiQg. From Eq. (2), 
the dominant part of 5^ is (6 ■ Vif / g, which in a chaotic flow grows exponentially with 
time. 

In Lagrangian coordinates, the induced current j — 'V x Bf/iQ is [3] 

„ a / b'' \ ,„, 



which after using Eq. (2) and projecting along the {u, m, s} basis becomes 

A*oju= A^6uU-VoXu-A^6„,s- [Volog(A^|6^|)-(m-Vo)m] +0(A^), (7) 
/^ojm= A^6„s- [Volog(A^|6u|)-(u-Vo)u]+0(Ay , (8) 
//ois = -Au bu m • [Vo log {Al\b^\) - (u • Vo)u] + O(A^) , (9) 

where ju := (s/gu-j), 6u := (u-6/y^), etc. The subscript on Vq denotes differentiation with 
respect to the Lagrangian coordinates, a. In Refs. [3, 4] it was concluded that the dominant 
term in g — A^ + A^ + A^j^ arises from A^ j^, so that overall the magnitude of the 
current grows as A^/g'^. This is an extremely rapid growth: following Ref. [3], we define the 
dissipation time scale, 

Td := Es/vf, (10) 

the ratio of the magnetic energy to the power dissipated in the plasma through Ohmic 
heating, and find that ~ 5'A~^. Of course, because the Cauchy solution for the magnetic 
field is only valid for ideal evolution — before the dissipation comes into play — the time 
scale Td must be regarded as an indicator of the accumulation of small scale gradients. 
When Td, which is initially large because of the small resistivity, becomes of the order of 
typical macroscopic time scales, the ideal evolution approximation becomes invalid. That Td 
evolves as gA~^ indicates that the ideal evolution must be abandoned rather quickly, and 
that relatively little magnetic energy has accumulated when this time is reached. We will 
see in the following section that the situation is not as hopeless as it appears: the actual 
domain of validity of ideal evolution is considerably greater than the calculation of this 
section indicates. 
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IV. THE EFFECT OF DIFFERENTIAL CONSTRAINTS 

The conclusions of Section III regarding the decay of Td was based on the assumption 
that u • Vo X u in (7) was a nonexponential function. In Refs. [6, 7] it was shown that this 
is not the case: in three dimensions we must have 

u- VoXu-A-2max(^,Ay — ^ 0, (11) 

as long as the flow is chaotic and there exists a contracting direction, s. Equation (11) is one 
example of a differential constraint in chaotic flows; it is realized with exponential accuracy 
in time. Because it will recur often, we define the factor 

/C:=max(^,A^); (12) 

For a steady incompressible flow, the growth or decay of fC is algebraic, not exponential, so 
in that case we can neglect /C compared to exponential factors to get an estimate of growth 
rates. 

The constraint (11) and others [sec below, Eqs. (13)-(14)] were derived in Ref. [6] using 
the condition that the Riemannian curvature of the metric tensor Qpq vanishes. In Ref. [7] the 
constraints were derived under more general assumptions, and they were shown to be present 
in arbitrary dimensions; their convergence rate was also obtained. The reason why (11) holds 
is not readily apparent; it is due to the overconstrained nature of the Lagrangian derivatives 
of the characteristic directions and coefficients of expansions [7]. All the constraints used in 
this paper have been verified numerically to a high degree of precision [5-7] for a variety of 
maps and flows. 

In three dimensions, there are two more constraints in addition to (11) [6, 7], 



u-(u- Vo)s + s- VologAu~ A, — ^0, (13) 

m • (m • Vo)s + s • Vo log A^ ~ max(As , A,VA^) — ' 0. (14) 

If we add (13) and (14), we find after some manipulation 

^ Vo • (x/5 s) - s • Vo log As - max(As , As/Am) — > 0; (15) 



This is a compressible {^/g ^ 1) version of the constraint used in Ref. [4], originally derived 
in two dimensions in Ref. [5]. In two dimensions it is the only independent constraint. 

In addition, if it happens that the ih direction (the middle eigenvalue) is a contracting 
direction {i.e., Am decreases exponentially for large times), then there is a fourth constraint, 

u-(u- Vo)m + m- VologAu~ Am — ^ 0, for A^ < 1. (16) 

Thus there are three constraints, (11), (13), and (14), if Am is a neutral or stretching di- 
rection, and four if Am is a contracting direction. For a generic three-dimensional chaotic 
flow satisfying the assumptions herein (that is, at least one stretching and one contracting 
direction), there can be no further differential constraints that involve only first deriva- 
tives of the characteristic directions and coefficients of expansion [7]. It is however possible 
that fiows with special symmetries admit additional relationships (algebraic, differential, or 
integro-differential) among the characteristic directions and coefficients of expansion. One 
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example is that for incompressible flows the coefficients of expansion obey the algebraic re- 
lationship Au AinAg = 1, as mentioned in Section II. We do not consider other relationships 
in the present work. 

Armed with the constraint (11), we find that the growth of the component ju, given 
by (7), is actually proportional to /C defined by Eq. (12), and not A^. This means that the 
dominant contribution to comes from (8), 

= AtAl [b^s- [Volog {Al\b^\) - (u ■ Vo)u])' + 0(A2 max(A^,/)) , (17) 

so that ~ AyA^jj/f?^ rather than the value A^/g^ derived in Section 111. With this new 
scaling the revised dissipation time scale becomes ~ gA~^A^ instead of fifA^ ^. Consider- 
ably more magnetic energy can thus be created before resistivity comes into play. In fact, 
for the typical case of Am nonexponential {i.e., for steady flows) the time for ideal evolution 
is roughly doubled. Of course because the production of magnetic energy is exponential in 
the ideal phase, this leads to an exponentially larger amount of energy. 

The modifications to the growth of the parallel current when the constraint is taken into 
account is even more drastic. By parallel current, we mean the part of j along the magnetic 
field, 

j| := U ■ B f/B' . (18) 

In Refs. [3, 4] it was concluded that j aligns with B because the dominant part of 
comes from j^, so that jj^ ~ A^/g'^. But if we apply the constraint (11), we find instead 
that jy ~ A^K? / g"^. This is a radically different growth rate: for a steady incompress- 
ible flow {i.e., Am nonexponential and g = 1), the ratio of the previous to the new result 
is A^/A^ = A^. The most important imphcation is that we now have j| <S j\, so that the 
current does not align with the magnetic field. In terms of length scales, we have that the 
scale of variation of B along itself, is much greater than the perpendicular variation, 
The dominant contribution to jy comes from a mixture of ju and jm, whereas the dominant 
contribution to comes from jm- 

Note that when we speak of alignment of two vectors we are expressing this in terms 
of their scalar product. The scalar product is invariant under coordinate transformations, 
including the transformation between Eulerian and Lagrangian coordinates. This means that 
the alignment occurs both in the Lagrangian and Eulerian frames. This is the power of the 
Lagrangian trajectory approach: the behavior of the vectors is easily derived in Lagrangian 
coordinates, and the results automatically apply to the Eulerian frame if expressed in terms 
of scalar quantities. 

V. A SINGULAR INITIAL CONDITION 

So far we have only made use of the constraint (11). In the present section we offer an 
application of the other constraints. If we apply the constraints (14) and (13) respectively 
to the ju and jm components of the current, equations (7) and (8), we find 

A*o Ju = A^ 6^ u • Vo X u - A^ 6m s • Vo log (Am 1 6m |) + A-^ O (Am , 1 , Ay ^) , (19) 
/^oim = Au&uS • Volog(A,|6,|) + ^ 0(Am, 1 , Ay , (20) 
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where we have made use of the fact that, asymptotically, the Lagrangian derivatives 



of log Au, log Am, and logg^ satisfy [7] 

s • Vo log Au = s • Vo log Au + 0(A, , AJA^ , AJA^) , (21) 

s • Vo log Am = s • Vo log Am + 0(As , As/ Am , Am/Au) , (22) 

s • Vo log^ = s • Vo log^+ 0(As , As/ Am) , (23) 



with Au, Am, and time-independent functions of a, unique up to a multiplicative constant. 
Equations (21)-(23) reflect the contraction of coordinates along the s direction, which leads 
to exponential convergence of Lagrangian derivatives; these convergence rates were derived 
in detail in Ref . [7] . 

From (20), we see that there are two special initial conditions for bu- The first is to simply 
let 6u = 0, but this is not relevant here because then the magnetic energy does not grow at 
the rate A^. A more interesting possibility is to let 

b-u^cr'K\ (24) 

where c is a constant and g and A^ are defined in (21) and (23). The g^^"^ factor appears 
in (24) because of the definition :— (u ■ h/ Then the first term in (20) vanishes, and 
we find that the leading-order behavior of the current is 

f^AlJ/g\ J:=max(K:^^A-2), (25) 

in contrast to A^A^/^i^ for a generic initial condition. For an asymptotic rate such as (25), 
the dissipation time scale (10) is ^ gj J] for a steady incompressible flow («.e.. Am nonex- 
ponential and gi = 1), we have ~ 1. Hence, for this special initial condition the dissipation 
time can be of order unity. We thus have an exponentially growing magnetic field whose 
growth can be sustained for a relatively long time before dissipation becomes important. 

As for the parallel current jj^, it remains unchanged for the singular initial condition and 
is equal to A\K? /g"^; its main contribution is from j^. Thus, asymptotically, the singular 
initial condition has the magnetic field and current aligned. (In terms of the scale of variation 
of the magnetic field, we have l\\ ~ ^±.) 

The existence of a singular initial condition is made possible by the constraints (11) 
and (13)-(14), but it is also crucial that the Lagrangian derivatives (21) and (23) converge 
to well-defined time-asymptotic values, which is the case for almost all initial conditions [7]. 
Otherwise it would not be feasible to ask for an initial condition of the form (24) since h is 
constant in time (until the diffusive regime is reached). 

We shall comment further on the singular initial condition in the Discussion section at 
the end of the paper. 



VI. HELICITY GENERATION 

In the present section we turn our attention to the production of magnetic helicity. The 
magnetic helicity V, :— A - B, where A is the magnetic potential with V x A = B, evolves 
according to 

^ + V • (t;7i) = -77 (2 j . B - V • (A X j)). (26) 
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(Some terms were absorbed by an appropriate choice of gauge.) Helicity can only be created 
or destroyed if 7^ 0. (We refer to the terms on the right-hand side of (26) as the hehcity 
generation terms, but they are not sign-definite and can also destroy helicity.) The second 
helicity generation term in (26) creates no net helicity, as it vanishes upon integration 
over the entire volume (for suitable boundary conditions). However, it can move helicity 
from one point to another: it represents a flux of helicity, so that the total helicity flux 
is vH — 7] A X j. 

The magnetic energy Eb = -B^/2/xo is dissipated by Ohmic heating in the plasma at a rate 
given by rjj'^. We want to examine the relative magnitude along fluid trajectories of local 
helicity generation as given by (26) and energy dissipation through Ohmic heating rjj^. This 
will help determine whether the net helicity in a region due to relative transport (helicity 
flux) or creation (the 2f]j-B term). It is thought that helicity generation is an important 
ingredient for the creation of a large-scale magnetic fleld [13], in part because helicity decays 
on a slower timcscale than energy [14] so its presence ensures a long-lived magnetic field. 
It is also important because of its implications for the topology (knottedness) of magnetic 
field lines [15-17]. Given the new results involving the constraints, is reasonable to ask if 
there is sufficient time in the ideal energy evolution phase to create a significant amount of 
helicity. Put another way, when the Ohmic heating term becomes of order one, how large 
are the helicity generation terms? 

Following the trajectory of a fiuid element, (26) becomes 



d_ 
di 



{A-b) = -r]{2j-b-Vo-{Axj)) (27) 



where all the variables are expressed in the Lagrangian frame and Ti = A ■ b/^/g. If we 
form the ratio of the helicity generation terms on the right-hand side of (27) and power 
dissipation r/j'^, the resistivity 1] cancels, making these two ratios convenient measures of 
the relative strengths of the effects. We can now ascertain whether the helicity generation 
terms become important before ideal evolution ends, and find which helicity generation term 
is largest. 

We start by comparing the first term on the right-hand side of (27) to the Ohmic heating 
term rjj^, ignoring the constraint (11): 

The helicity generation term is typically exponentially small compared to energy dissipation. 
Including the constraint (11) does not improve matters: 

so that the helicity generation term is still typically exponentially smaller than Ohmic heat- 
ing. 

For the singular initial condition of Section V, the ratio of helicity generation to Ohmic 
dissipation yields 

j ■ b KKi K 9^ gJC 
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The ratio is nonexponential {i.e., order one) for a steady incompressible flow. This opens 
the possibihty of creating a sizable amount of helicity before Ohmic dissipation becomes 
important. 

The second term on the right of (27) — the flux of magnetic helicity in Lagrangian 
coordinates — is given by 

3 

Vo • ( A X j) = 5] Vo • [Al (A X e^) jj^] (31) 

where // = {u, m, s} and {e} := {u, ih, s}. Because A does not acquire a special orienta- 
tion with respect to e^, upon expansion the divergence in (31) will generate terms of the 
form u • Vo log A, and these grow proportionally to Au [7] . Hence we flnd 

3 

Vo-(Axi)~^^Ajj,. (32) 

The overall growth rate thus depends on whether the constraint (11) is used, because this 
affects which component of is largest. Without the constraint we have both ju and 
proportional to A^/ y^; then the dominant contribution in (32) comes from j^, yielding 

Vo-iAxj)r.^Alj^r^Al/g. (33) 

But the result (33) was obtained without taking differential constraints into account and is 
thus incorrect. Upon applying the constraint (11), we have ~ ^/s/g, so that 

Vo-{Axj)r^Al)C/g. (34) 

Taking the ratio with the Ohmic heating 77 j^, we flnd that Vq ■ {A x is proportional 

to g without the constraint (11), and to gA~^A'^ /C with the constraint. For the case of 
a steady incompressible flow, this helicity source term is again exponentially smaller than 
the power dissipation, and moreover the ratio of the two is roughly A~^ both with and 
without the constraint. The rate (34) also applies to the singular initial condition, but in 
that case j'^ grows more slowly [Eq. (25)], allowing more accumulation of helicity before the 
end of ideal evolution. 

If we compare the two terms (helicity creation and helicity flux) on the right-hand side 
of (27) to each other, we flnd 

^^^4^ ~ A„ (36) 

in all three cases (constrained, unconstrained, and singular initial condition). We conclude 
that though the proper application of the constraints changes the absolute magnitude of 
each hehcity source term in (27), when expressed as a ratio with the current their relative 
magnitude is unchanged. The helicity flux (second term) always tends to dominate the 
helicity creation (flrst term) as a helicity source on the right-hand side of (27). 

Table I summarizes the growth rates derived in previous sections and in this section 
for helicity generation. For a clearer picture of relative magnitudes. Table II gives the 
corresponding growth rates for a steady incompressible flow. The singular initial condition 
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No constraint Constraint Sing, initial cond. Comments 





K/9 


K/9 


K/9 


Magnetic energy 


f 




KJ^l/9' 


AIJ/9' 


Power dissipation 


•^11 


K/9' 


Alley 




Parallel current 


j ■ b/f 


9K' 


2 A-2 /c 


gic/J 


Helicity creation vs power 


Vo • (A X 


9K' 




gkJZlJ 


Helicity flux vs power 



TABLE I: Asymptotic growth of physical quantities. Three cases are considered: (i) No differential 
constraint; (ii) Constraint given by Eq. (11) used; (iii) Singular initial condition of Section V, which 
uses the constraints (11)-(14). The factor /C is defined in (12), J in (25). 





No constraint 


Constraint 


Sing, initial cond. 






AS 


A2 




A^ 


A^ 


A2 




A6 


h? 


A2 


3 ■ bif 






1 


Vo • (A X 






Au 



TABLE II: As for Table I but assuming a steady incompressible flow (Am ~ 1, g = 1). 

is the only case that allows for the possibility of significant helicity creation or flow before 
the dissipative regime is reached. For the other two cases (unconstrained and constrained), 
for very small rj the magnetic field will have built up huge gradients by the time the evolution 
ceases to be ideal. Exponentially large amounts of power thus seem required to create or 
move helicity. Even for the singular initial condition, there is no guarantee that the helicity 
will be created rather than destroyed, or will be concentrated rather than dispersed, since 
the generation and flow terms are not deflnite in sign. In fact, it is known that substantial 
cancellation of the helicity occurs [18], a fact that is not captured in our analysis. 

We emphasize that the results for the unconstrained case are presented merely for com- 
parison purposes: the constrained results always apply because differential constraints are 
an unavoidable geometrical consequence of the chaotic nature of the flow [7] . The singular 
initial condition, however, only applies to situations that are initialized in the special manner 
given by (24). Otherwise the constrained result apphes. 

VII. DISCUSSION 

Using Lagrangian coordinates, we have revised earlier estimates of the magnitude of 
energy dissipation (Section IV) and hehcity generation (Section VI). The estimates were 
corrected for differential constraints on the rates and directions of stretching. The inclusion 
of constraints leads roughly to a doubling of the period of ideal evolution, for the case of a 
steady incompressible flow. We also found after using constraints that he current does not 
align with the magnetic held in the ideal phase, in contrast with the results of Refs. [3, 4] 
where it was concluded that it did. The constraints thus have a deep implication on the 
evolution of the magnetic field, at least in the ideal phase. 

In spite of the constraints, we find that the helicity generation terms are exponentially 
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smaller than the energy dissipation, so that creation of significant amounts of hclicity cannot 
occur in the period of ideal energy growth (when resistivity can be neglected). Our analysis 
does not take into account the fact that the helicity created will be fractal in nature, and so 
very little of it will contribute to a large-scale magnetic field, as pointed out by Gilbert [18] . 
This suggests a close inspection of the energy dissipation in numerical simulations, since 
current dynamo models could require vast amounts of energy to create helicity. 

We found (Section V) that there exists a special initial condition with desirable proper- 
ties: (i) the power dissipated is of the same order as the magnetic energy; (ii) the hehcity 
generation terms are comparable to the power dissipated. This opens the possibility of a 
fiow creating a considerable amount of energy and helicity before reaching the end of the 
ideal evolution regime. The physical relevance of the singular initial condition is uncertain. 
We cannot expect a realistic system to have an initial magnetic field corresponding that 
particular condition. The part of the initial condition that is not singular will reach the dis- 
sipation scale much sooner, so the ideal evolution hypothesis becomes invalid. But does the 
singular part of the initial condition survive, or is it dissipated away? An analysis involving 
the full, dissipative equation is necessary to establish this and has not yet been performed. 
If it does survive, it could be closely related to the "strange eigenfunction" (i.e., a fractal) 
that occurs in the freely-decaying advection-diffusion problem [19] and its analogue for the 
dynamo [1]. 
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